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We explore the potential of double core hole electron spectroscopy for chemical anal- 
ysis in terms of x-ray two-photon photoelectron spectroscopy (XTPPS). The cre- 
ation of deep single and double core vacancies induces significant reorganization of 
valence electrons. The corresponding relaxation energies and the interatomic relax- 
ation energies are evaluated by CASSCF calculations. We propose a method how 
to experimentally extract these quantities by the measurement of single and double 
core- hole ionization potentials (IPs and DIPs). The influence of the chemical envi- 
ronment on these DIPs is also discussed for states with two holes at the same atomic 
site and states with two holes at two different atomic sites. Electron density differ- 
ence between the ground and double core-hole states clearly shows the relaxations 
accompanying the double core-hole ionization. The effect is also compared with the 
sensitivity of single core hole ionization potentials (IPs) arising in single core hole 
electron spectroscopy. We have demonstrated the method for a representative set 
of small molecules LiF, BeO, BF, CO, N 2 , C 2 H 2 , C 2 H 4 , C 2 H 6 , C0 2 and N 2 0. The 
scalar relativistic effect on IPs and on DIPs are briefly addressed. 
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I. INTRODUCTION 



The effect of the chemical environment manifests itself in energy differences of molecular 
core levels with respect to the atomic ones referred to as "chemical shifts". These can be 
measured by core level spectroscopies, e.g., by x-ray photoelectron spectroscopy (XPS) also 
known as electron spectroscopy for chemical analysis (ESCA) and by x-ray-induced Auger 
electron spectroscopy (XAES) 1 . Both spectroscopies have shown to be exceedingly successful 
tools to reveal the quantitative elemental composition of molecules and solids. 

More than two decades ago, Cederbaum et al.- 1 ^ discovered that the creation of double 
core vacancies in molecular systems probes the chemical environment more sensitively than 
the creation of single core vacancies. Two-atomic site double ionization potentials, or briefly 
two-site DIPs ( or two-site double ionization energies, DIEs ) are particularly sensitive to 
the chemical environment as the examples of the C2H2, C2H4, C2Hg^ and CqB.^- molecules 
demonstrate. The chemical shifts of one-atomic site DIPs, or briefly one-site DIPs, were 
found to be similar to the chemical shifts of the single core level ionization potentials (IPs), 
or ionization energies (IEs). This finding has given impetus to a number of theoretical 
studies aimed at elucidating properties of molecular double core hole states^-. 

So far experimental explorations of double core hole states with conventional XPS were 
restricted to those having two vacancies on the same atomic site only^ since the probability 
to produce a two-site double core hole state with one-photon absorption is practically zero at 
third-generation synchrotrons due to low x-ray intensities. This prevented further progress 
of the subject. The situation has changed with development of x-ray free-electron lasers 
(x-ray FELs) 10 . At FEL facilities in operations, such as FLASH in Hamburg 11 and SPring-8 
Compact SASE Source (SCSS) test accelerator—, multi-photon absorption processes result- 
ing in multiply ionized states of various systems have been extensively studied^r—. In the 
x-ray FEL facility LCLS at SLAC National Accelerator Laboratory, which has just started 
operations^, ultrashort pulses of a duration about 1-fs containing 2.4xlO n photons with 
energies of 1 keV are expected to be generated^>22 thus opening up the possibility to study 
molecular two-site double core hole states. Inspired by the advent of the x-ray FEL at 
LCLS, Santra et al.— have demonstrated theoretically by the proof-of-principle simulations 
on the organic para-aminophenol molecule that two-site double core hole states can indeed 
be probed by means of x-ray two-photon photoelectron spectroscopy (XTPPS). 
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The operating principle of XTPPS is depicted schematically in Fig. [TJ The initial step 
in XTPPS corresponds conventional XPS, i.e., a neutral molecule with an energy Eq is 
irradiated by an x-ray photon with an energy wj and a photoelectron with the kinetic 
energy k pl /2 is ejected. This photoelectron carries information about a singly core ionized 
state E + of the molecule. If a second x-ray photon is absorbed before the intermediate core 
hole state decays, the second photoelectron expelled from the cation with the kinetic energy 
kp 2 /2 carries information about a double ionization potential. It is important to have an 
intense x-ray pulse with a duration that is significantly shorter than the core-hole lifetimes 
(typical lifetimes of core ionized states of F, O, N and C atoms are 3 to 7 femtoseconds). 
If the pulse duration is longer than these lifetime, then Auger decay is likely to occur prior 
to absorbing the second photon and thus the double core hole states may not be probed. A 
dicationic state E ++ of the system prepared by two-photon absorption decays electronically. 
Two primary Auger decays take place which overlap in time. An Auger decay happens 
preferably at that atomic site where the core hole has the shorter lifetime and an Auger 
electron with kinetic energy k\ 1 /2 is ejected. This process proceeds in the presence of the 
second core hole which also decays via the Auger mechanism emitting an electron with kinetic 
energy k^ 2 /2. The electrons ejected via such a cascade of Auger decays can in principle 
be measured by a novel Auger spectroscopy which we call x-ray two-photon-induced Auger 
electron spectroscopy (XTPAES). 

It is worthwhile to note that double core ionization can be accompanied by various shake- 
up processes similar to single core ionization. These many-body effects should manifest 
themselves in XTPPS spectra as satellites which are of interest as well. Both x-ray two- 
photon-induced Auger spectra and satellites structures will be addressed elsewhere. 

The subject of the present paper is the main double core hole states. In order to provide 
a guideline for XTPPS experiments, we have performed ab initio calculations of core level 
single and double ionization potentials of LiF, BeO, BF, CO, N 2 , C2H2, C2H4, C2H 6 , CO2, 
and N 2 molecules. In addition we have explored the sensitivity of the DIPs to the chemical 
environment of the core ionized atoms. We decompose the DIPs in three physical contribu- 
tions such as the orbital energies, the electrostatic repulsion energy between two core holes 
and the generalized relaxation energy and describe how the latter can be extracted from the 
experimental XTPPS spectra. 
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II. COMPUTATIONAL METHODS 



Ab initio calculations of the vertical ionization potentials of the single and double core 
vacancy states of LiF, BeO, BF, CO, N 2 , C2H 2 , C2H4, C2H 6 , CO2, and N 2 were per- 
formed using the ASCF— and CASSCF 23 methods. The molecular geometries used in these 
calculations were optimized at the M0ller-Plesset level of theory (MP2)^ employing the 
correlation-consistent polarized valence triple zeta (cc-pVTZ) basis sets of Dunning^. De- 
pending on whether singly or doubly ionized states were considered, the configurations in the 
CASSCF method were restricted to those having one or two holes in the K-shell orbitals, re- 
spectively. We used the active spaces comprising all the occupied molecular orbitals (except 
for the Is orbitals of the atoms other than H) and all the valence unoccupied n* and a* ones 
which contain large contributions from different atomic p orbitals. Thus, the active space of 
the CASSCF calculations consists of 2s, a, tt, it* and a* orbitals with core occupancy being 
fixed. The cc-pVTZ basis sets were employed in all the CASSCF and ASCF calculations. 
For CO and C 2 H 2 , the cc-pVDZ, cc-pCVTZ and cc-pVQZ basis sets were also used in order 
to examine the basis set dependence of our results. All the single and double core-hole states 
were solved by independent CASSCF calculations using different configuration space and, 
therefore, the calculated states are not strictly orthogonal to each other. However, the one 
and two-site double core-hole states are well separated in energy and their interaction are 
expected to be negligible. 

For molecules with equivalent atoms, N 2 , C2H2, C2H4, C2H 6 and CO2, we calculated 
double ionization potentials using both localized and delocalized molecular orbital pictures 
following the recipe given by Cederbaum et al.— and discuss differences between them. Note 
that only core orbitals were localized which was performed with the Boys method^. In the 
localized representation we obtained the ionization potentials of the one-site double core hole 
states S^~ 2 and S^ 2 as well as of the two-site double core hole states Sj" 1 S 2 ~ 1 . Carrying out 
calculations with wave functions described by the linear combinations S|f 2 ±S^ 2 gives rise to 
double ionization potentials in the delocalized picture. Differences between single ionization 
potentials arising due to applying localized and delocalized representations are not studied 
in the present paper because they have been discussed in detail before 2 ^ 1 ^. 

In this work, we ignore the geometry relaxation of ionized state. In core ionization 
the change of geometry can be significant, depending on the case under investigation. In 
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single core ionization one can explain the measurements well by employing the concept of 
vertical transitions. As in single ionization, also in XTPPS where the two X-ray photons 
must be absorbed within a time shorter than the Auger decay times, the concept of vertical 
transitions can be expected to be very useful. 

In order to assess the impact of scalar relativistic effects on the core level single and double 
ionization potentials we made relativistic CASSCF calculations for the CO and BF molecules 
using the eighth order Douglas-Kroll-Hess Hamiltonian (DKH8) 29 ~ 33 . To get insight into the 
dynamic correlations, we also performed CI calculations with the CAS space plus single 
excitations from the CAS for both single and double core-hole states. 

All calculations were done with the Molpro2008 quantum chemistry package 3 ^. 



III. RESULTS AND DISCUSSIONS 

A. Single core hole states 

Let us first discuss single core hole IPs. The ionization potential for the formation of a 
vacancy S^ 1 can be represented as 

IP(S- 1 ) = -e s -RC(S- 1 ), (1) 

where Es is the corresponding orbital energy and i?C(S'~ 1 ) is a contribution to the ionization 
potential due to relaxation R^S' 1 ) and correlation C(5' _1 ) effects: 

RCiS- 1 ) = R^S- 1 ) + C^S- 1 ) (2) 

The relaxation and electron correlations intermix with each other and cannot be strictly 
separated. The separation of these quantities was discussed in details in a perturbative 
way2£ and in a nonperturbative way using MRCC 36 . The correlation contribution can be 



37 and 
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where CI describes a 



further decomposed into two parts CI and C2 (see Refs. 
part of the ground state pair correlation energy disappearing upon removal of an electron 
from the spin orbital S, and C2 accounts for changes in the remaining pair correlation 
energy due to relaxation. Except for CI which is a very small contribution, all contributions 
to RC are thus associated with relaxation of molecular orbitals. Therefore, for brevity of 
discussion, we may call RC the generalized relaxation energy. 



A straitforward way to obtain the relaxation energy is to perform ASCF calculations. 
i^S 1 " 1 ) is then derived as the difference between the respective orbital energy (with opposite 
sign) and the calculated IP. In order to get a correlation contribution to IP, post-Hartree- 
Fock calculations are generally needed. CASSCF is one of these methods. Noteworthy, in 
systems with core holes delocalized due to symmetry requirements, C2 can be accounted by 
performing ASCF calculations using localized orbitals instead of delocalized ones as shown 
by Cederbaum and Domcke^. 

In Table [I] we list IPs calculated with the CASSCF and ASCF methods together with 
available experimental values^— . Table [T]also contains the constituting parts of IPs, namely 
the orbital energies, the relaxation energies obtained from ASCF calculations, as well as 
the generalized relaxation and pure correlation contributions, both obtained from CASSCF 
calculations. The correlation contributions were calculated by subtracting the CASSCF IPs 
from the ASCF ones. Note that, since ASCF calculations for molecules with equivalent 
atoms were performed using the localized representation, the calculated IPs correspond to 
the localized Is orbitals rather than to the delocalized lo~ g and la u . 

In general, the agreement between the CASSCF and experimental results is reasonable. 
This concerns both the absolute values of IPs and the g-u energy splittings for molecules 
with equivalent atoms. Except for basis set effects which always are an issue in ab initio 
calculations, and relativistic effects, deviations from the experiment are attributed to the lack 
of dynamic correlations in the ground and single core hole states, and to the core-valence 
separation approximation employed in the calculations. We notice that influences of the 
above-mentioned effects and approximations partially compensate for each other. Indeed, 
performing calculations without core- valence separability lowers IPs 45 . A lowering of IPs can 
also be achieved by improving basis sets. On the other hand, taking into account relativistic 
effects increases IPs. In Appendices |A] and |B] we explore the basis sets and relativistic effects 
in more detail. 

It is interesting to compare the different contributions to IPs in Eq. (JT|). After the 
orbital energy, the relaxation energy represents the largest constituent part of a core level 
IP. It increases nearly proportional to the atomic number Z. For some molecules, however, 
remarkable deviations from this trend occur under influence of the chemical environment. 
A crucial role for the relaxation energy plays the change of the electron density distribution 
of valence electrons, Ap, in an atom due to a formation of chemical bonds with neighbors, 
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and the interaction of Ap with the core hole. Ionic bonds give rise to the strongest changes 
of the electron density distribution. As a consequence, the relaxation energies associated 
with core ionization of electron acceptors in ionic molecules (e.g. O and F in BeO and LiF, 
respectively) are noticeably larger than the relaxation energies of the same atoms bound 
by covalent bonds with their neighbors (O and F in CO and BF, respectively). For other 
factors influencing the relaxation energies see Ref. |2|. 

In comparison to relaxation effects, correlation effects induced by core ionization are 
rather small. According to our calculations, the magnitude of the static correlation effects 
does not exceed 3 eV for the molecules studied and accounting for missing dynamic corre- 
lation can hardly modify this situation dramatically. Interestingly, the largest correlation 
effects manifest themselves in atoms whose neighbors are the strong electron acceptors O 
and F. 

The effect of the chemical environment on core level ionization potentials of various 
systems including the molecules explored here is rather well established and we therefore 
refrain from long discussions on this subject. We only mention that the chemical environment 
is able to introduce large changes in the ionization potentials as, for example, can be realized 
by comparing molecules with ionic and covalent bonds. On the other hand, in the sequence 
of the C2H2, C2H4, C2H6 molecules characterized by the triple, double and single carbon- 
carbon bond, respectively, the impact of the chemical environment is rather moderate. In 
contrast to single core hole ionization potentials, double core hole ionization potentials reveal 
much more pronounced sensitivity to the chemical environment as it was first demonstrated 
in Refs. Sandfl 

B. Double core hole states 

1. General equations and results 

In analogy to Eq. ([[]), we represent the double ionization potential of a state with two 
core vacancies S 1 , -1 and S~ l as 

DIP{S~\ S7 1 ) = -e Sl - e Sl - RC(Sr\ S^) + E RE , (3) 

where Ere is the repulsion-exchange energy of the two core holes. For an one-site double 
core hole state, it is described by the two-electron integral Vs i s,s i s t : or (SiSi\SiSi), and, 
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for a two-site double core hole state, by a linear combination of the integrals Vs^s^ and 
VsfSjSjSi where the exchange term is negligibly small when the core holes are well localized 5 . 
The generalized relaxation RC(S^ 1 , S^ 1 ) can be decomposed into three parts 

RC{S'\ Sj 1 ) = RC(S; 1 ) + RCiSj 1 ) + NRC(S-\ S" 1 ), (4) 

where RC^Sj 1 ) given by Eq. (TS]) describes relaxation and correlation effects induced by 
creation of the core vacancy Sj~ as there were no core vacancy Sj . The relaxation and 
correlation energies are expected to be non-additive upon creation of multiple vacancies. A 
possible deviation from additivity is thus described in Eq. 03]) by the non-additivity term 
NRCiS- 1 ^; 1 ). 

Depending on whether two core holes were created on the same atomic site or on dif- 
ferent atomic sites, NRC(S~ 1 , Sj ) may be called the excess generalized relaxation energy, 
ERC^Sj 1 , Sj -1 ), or the interatomic generalized relaxation energy, I RC '(Sj -1 , Sj 1 ) . Note 
that, while the RC{Sj x ) and ERC{Sj\ Sj 1 ) measure local properties of a core ionized 
atom, IRC{Sj 1 ,Sj 1 ) measures the impact of the environment "between" the atoms in- 
volved. 

In Table [TT] we list the calculated double ionization potentials. We also show the correla- 
tion contributions to DIPs. As one can see these contributions are remarkably larger than 
those to the single IPs and may constitute 5.6 eV. In the special cases, however, when we per- 
formed calculations with delocalized core orbitals, differences between ASCF and CASSCF 
values rise to 27-35 eV resulting from the failure of the ASCF method in the delocalized 
picture to account for all relaxation contributions as described by Cederbaum et al.->2. 

One can notice by comparing Table [land [IT] that the impact of the chemical environment 
is different for double and single ionization potentials. Of particular interest is to compare 
two-site DIPs with single IPs since their sensitivities to the chemical environment reveal 
major differences. A prominent example already discussed in detail in Ref. 2] is the hydro- 
carbons C2H 2 , C2H4 and C2H 6 . Here, the chemical shifts in the two-site DIPs are much 
more pronounced than in the single IPs being also attributed to different carbon-carbon 
bondlengths resulted from a different number of hydrogen atoms in these molecules. In XPS 
one can hardly distinguish between these three compounds while in XTPPS this should 
be possible in principle. The situation is somewhat different for the individual molecule 
N 2 0. In this molecule the sensitivity of two-site DIPs to the chemical environment is lower 



compared to that of single IPs. Indeed, the core ionization potentials of the terminal and 
central nitrogen atoms differ by 4 eV whereas the difference between the Ntls _1 01s _1 and 
N c ls _1 01s _1 double ionization potentials constitutes 2.3-2.8 eV. The latter energy difference 
is much lower than 11 eV which one would expect taking into account only the differences 
between the N c O and N t O bondlengths and between the single core hole ionization poten- 
tials. As we show below, the reason for such a dramatic reduction of the chemical shift has 
to do with distinct relaxation processes induced by the creation of different pairs of core 
holes. 

Taking into account Eqs. ([I]) and (@J, we can represent DIP(Sj x , Sj ) as 

DIP(Sj\ Sj 1 ) = IPiSr 1 ) + IPiSj 1 ) - NRC(Sj\ Sj 1 ) + E RE , (5) 

and define the ionization potential of the core vacancy Sj in the presence of the core 
vacancy Sj 1 as 

IP{Sj\ Sj 1 ) = IP^Sj 1 ) - NRC(Sj\ Sj 1 ) + E RE , (6) 

whereas 

IP(Sr l ; Sj 1 ) = IP(Sr l ) - NRC(Sj\ Sj 1 ) + E RE (7) 

is defined as the ionization potential of the core vacancy Sj in the presence of the core 
vacancy Sj 1 . 

Both the ionization potentials IP^Sj 1 ) of a neutral system and the ionization poten- 
tials IP^Sj 1 ; Sj 1 ) of a core- ionized one can be obtained experimentally, e.g., by means of 
XTPPS. In XTPPS, the kinetic energy KE^Sj 1 ) of the first photoelectron ejected from the 
orbital Si defines IP^Sj 1 ), whereas the kinetic energy KE^Sj 1 ; Sj 1 ) of the second photo- 
electron ejected from the orbital Sj defines IP(Sj ; Sj 1 ). Obviously, the sum IP^Sj 1 ) + 
" 1 ) gives DIP(S i ,Sj 1 ). As shown below, important properties of the system 
under study can be extracted also from the measurable energy difference 

AE = KE(S^) - KE{SJ\ Sj 1 ) = IP{Sj\ Sj 1 ) - IP^Sj 1 ) 

= DIP(Sj\ Sj 1 ) - IPiS- 1 ) - IPiSr 1 ). (8) 

Similar to DIP(Sj 1 , Sj 1 ), the kinetic energy KE(Sj 1 ; Sj 1 ) depends significantly on 
the mutual arrangement of the core vacancies Sj and Sj 1 in a molecule. This is clearly 
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seen from Table II III where we collect the kinetic energies of all the core electrons of the 
CO molecule which one would detect in an XTPPS experiment given that the molecule 
is irradiated by an x-ray pulse with photon energies of 1 keV. First of all, we notice that 
it is more difficult to remove an electron from the core ionized CO molecule than from 
the neutral one. The respective energy difference is about 70-90 eV when the first and 
the second core electrons are ejected from the same core orbital. This energy difference 
reduces drastically to about 15 eV when different core orbitals are affected. Apparently, the 
electrostatic interaction between the two core holes plays a crucial role here. NRC(S~ 1 , S^ 1 ) 
exerts an influence on the above energy differences too, as can be deduced from Eqs. (EJ and 
©■ 

2. One-site double core hole states 

If Si = Sj = S then AE takes the form 

AE1(S' 2 ) = -ERC(S-\ S- 1 ) + V ssss . (9) 

We calculated AE1(S~ 2 ) for the molecules under study using the respective CASSCF single 
and double core hole ionization potentials and collect them in Table IIVI The dependence of 
AE1(S~ 2 ) on the atomic number Z is displayed in Fig. |2Ja). 

The excess generalized relaxation energy ERC{S~ 1 , S 1-1 ) can be easily obtained by mea- 
suring the energy difference AE1(S~ 2 ) provided that the integral Vssss is known. Vssss 
can be extracted from ab initio Hartree-Fock calculations on the electronic ground state of 
neutral molecules. Alternatively, it can be calculated by using the approximate analytical 
expression suggested in Ref. [2] : 

Vssss = (2 5 / 2 /3tt)(Z-2- 3 / 2 ). (10) 

The respective results for Vssss an d a difference between them are discussed in Appendix [Cl 
Using Eq. @ we represent ERC^S^ 1 , S^ 1 ) as 

ERC(S-\ S- 1 ) = RC(S-\ S- 1 ) - 2RC(S- 1 ). (11) 

It has been shown in Ref. h| that at the second order perturbation theory the following 
relationship between the relaxation energies is valid: 

R(S- 1 ,S- 1 )=AR(S~ 1 ). (12) 
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Since the impact of correlation into ionization potentials is small compared to the impact of 
relaxation, we expect that a similar relationship exists between the generalized relaxation 
energies RC(S~ 1 , S^ 1 ) and RC^S^ 1 ). Let us therefore introduce that 

RCiS- 1 ^- 1 ) = nx RCiS- 1 ), (13) 

and find the optimal n. After the substitution of Eqs. ffTTj) and f fT3"j) into we get 

AE1(S- 2 ) = -(n - 2)RC(S- 1 ) + V SS ss- (14) 

Now we can easily calculate n by using the ab initio results for AE1, RC and Vssss- The 
respective values of n as a function of the atomic number Z are shown in Fig. [3j As one 
can see, deviations of the calculated n from the expected value of 4 are rather small (15% 
in the worst case of Li) and therefore n = 4 can be considered as a plausible approximation 
for the molecules studied in the present work. 

As a result, we obtain the following expression for the generalized relaxation energy: 

RCiS- 1 ) = (V S sss ~ AEl(S- 2 ))/2. (15) 

The values of RC^S' 1 ) calculated by means of Eq. ( Tl5|) are given in Table [TV] and also 
plotted in Fig. Mjo) as a function of Z. It is worthwhile to note a reasonable agreement 
between them and the corresponding results from Table [H 

3. Two-site double core hole states 

If Si Sj then AE takes the form 

AE2(S-\ S- 1 ) = -IRC(S-\ S- 1 ) + 1/r, (16) 

where the repulsion-exchange energy Ere has been approximated by the inverse of the 
distance r between the two atoms with the core vacancies S^ 1 and S~ x . We calculated 
AE2(S^ 1 , Sj) using the CASSCF double and single core hole ionization potentials and 
collected these results in Table [TV] 

By looking at Eq. ffl6|) . one expects a linear dependence between AE2(S^ X , Sj~ ) and 
1/r . This expectation is not realized however as seen from Fig. IH(a) where a variation 
of AE2 with 1/r is shown. To elucidate the complex behavior of AE2 we calculated the 
interatomic generalized relaxation energy IRCiS' 1 ^' 1 ) from Eq. (JE) and plotted these 
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results as a function of r in Fig. BJb). We found both positive and negative values of IRC 
(see also Table IIVI) which indicate on an enhancement or suppression of relaxation effects, 
respectively. 

In the case of diatomic molecules, IRC is always negative and thus the relaxation is 
suppressed. This suppression may be interpreted by considering the change of the electron 
density. A core hole S~ l attracts valence electrons and increases the electron density in its 
vicinity, yielding a deficiency of the electron density in the vicinity of the atom with a core 
orbital Sj. The relaxation energy corresponding to the creation of the core vacancy S 1 " 1 in 
the presence of the core hole S"" 1 is therefore reduced. The amplitudes of IRC are smaller for 
LiF, BF, and N 2 than for BeO, and CO. This is because the electrons are strongly attracted 
to the F site in LiF and BF, or tightened in the triple bond in N 2 and thus the electron 
density flow due to core hole creation is suppressed by these chemical environments. 

Values of IRC(S^ 1 , S' 1 ) for the polyatomic molecules C2H2, C2H4, and C2H 6 are, in 
contrast, positive. The enhancement of the relaxation for these molecules occurs due to the 
electron density flowing from the C-H bonds to the two carbon core hole sites. 

The interatomic generalized relaxation energy exhibits a very interesting behavior with r 
in triatomic molecules as seen from Fig. @|b). IRC(S^ X , S" 1 ) is positive for C0 2 and N 2 
when S^ 1 and S^ 1 are in adjacent atoms, namely in C and O in C0 2 and in N t and N c or in 
N c and O in N 4 N c O. In these cases, the third atom plays the role of an electron donor and 
enhances the relaxation of the double core hole in the other two sites. The values of IRC 
are, on the other hand, negative for C0 2 and N 2 with two holes in the terminal atoms, 
namely with holes in the two O sites of C0 2 and with holes in the N t and O sites of N 2 0. In 
these cases, the creation of the core hole on one site already withdraws the electron density 
from the central atom and thus reduces the possibility of relaxation due to the creation of 
the second hole in the other terminal site. 

In order to analyze the reorganization caused by double core hole ionization, we calculated 
the electron density difference between the ground and double core hole ionized states. 
These electron density differences without the Is contribution are plotted in Fig. [51 to 
better visualize the reorganization of the valence electrons. In the blue or green region, the 
electron density increases, while the density decreases in the red region. In the Cils~ 2 state 
of C 2 H 4 , the electron density of the Ci-C 2 , Ci-H bond and H atoms connected to Ci atom 
reduces and flows to the region around Cj atom as shown in Fig. Ufa). In the Cils _1 C 2 ls _1 
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state of C2H4, on the other hand, the electron density in the region of C-H bonds and H 
atoms flows to the region of both C atoms as in Fig. [5](b) . This explains the positive value 
of IRC in the Cils^C^ls -1 state as noted above. In the case of the Ols" 2 state of CO, 
the electron density of the CO bond is used for the reorganization around the O atom (Fig. 



IV. CONCLUSIONS 

We have computed the ionization potentials of single and double core hole states of the 
small molecules LiF, BeO, BF, CO, N 2 , C 2 H 2 , C 2 H 4 , C 2 H 6 , C0 2 , and N 2 by means of 
the ASCF and CASSCF methods in order to explore the impact of the chemical environ- 
ment on the respective ionization processes and provide a guidance for x-ray two-photon 
photoelectron spectroscopy (XTPPS) experiments. 

Our calculations have demonstrated that except for N 2 0, the double ionization poten- 
tials, especially the two-site ones of all these molecules are more sensitive to the chemical 
environment than the single ionization potentials. The sensitivity to the chemical environ- 
ment of the two-site DIPs of N 2 is mainly governed by the interatomic relaxation energies 
which, in turn, strongly depend on the electron density distribution between the core-ionized 
atoms. 

The quantities extracted from XTPPS are differences between the kinetic energies of core 
electrons ejected via the first and second ionization steps, i.e., of core electrons ejected from 
neutral and core-ionized systems, respectively. These kinetic energy differences are defined 
by a localization of the two core vacancies created and by relaxation processes induced 
by double core ionization. We have shown how one can extract the generalized relaxation 
energy associated with single core ionization as well as the excess and interatomic generalized 
relaxation energies associated with one-site and two-site double core ionizations, respectively, 
from experimental data by knowing the repulsion energy between the two core holes. The 
corresponding XTPPS experiments are now in preparation in the x-ray free electron laser 
facility LCLS at SLAC National Accelerator Laboratory. 
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Appendix A: Scalar relativistic effects 

It is essential in view of their large impact to take into account relativistic effects when 
considering systems with heavy atoms. For light atoms, relativistic effects are of less impor- 
tance but still should be accounted when highly accurate results are needed. In the present 
work we assess the impact of scalar relativistic effects on IPs and DIPs by carrying out 
calculations with the relativistic DKH8 Hamiltonian for the CO and BF molecules. The 
results of these calculations are shown in Table |VJ As one can see the magnitude of the 
scalar relativistic effects on the single ionization potentials grows with the increasing atomic 
mass constituting 0.03, 0.09, 0.35, and 0.59 eV in the case of the B, C, O and F atoms, 
respectively. One of us has shown that a similar tendency is observed also for the third-row 
Si, P, S, CI atoms^. A growth of scalar relativistic effects with the atomic mass exhibits also 
in the case of double core hole ionization. We note that the relativistic effects on the one-site 
DIPs are about 2.3-3 times larger than those on the respective single IPs. Interestingly that 
in the case of two-site doubly ionized states the relativistic effects are perfectly described by 
the sum of the relativistic effects associated with the constituting single core vacancies. 

Appendix B: Basis set effects 

In this section we explore the basis set dependence of the single and double core hole 
ionization potentials by examples of the C2H 2 and CO molecules. We have examined four 
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basis sets: the Dunning's correlation-consistent basis sets cc-pVXZ (X=D,T,Q) and the 
cc-pCVTZ one. The latter contains tight basis functions which are added for a better 
description of properties of core-level states. The results are collected in Table IVIl By 
looking at the sequence of the cc-pVXZ results we find significant differences between the 
cc-pVDZ and cc-pVQZ values, especially in the case of the one-site DIPs where differences 
of nearly 10 eV are present. The major changes occur upon improving the basis set from a 
double- to triple-zeta quality. Choosing the cc-pVQZ basis sets lowers ionization potentials 
by only 0.2-1.2 eV relative to the cc-pVTZ values. The cc-pCVTZ results are lower in energy 
than the cc-pVTZ ones whereby they nearly coincide with the cc-pVQZ results in the case 
of the CO molecule. We used the cc-pVTZ basis throughout as a compromise between the 
accuracy of the results and the computational costs. For the double core-hole states, the 
relaxation of valence orbitals is important, in particular for the one-site states where all Is 
electrons on one atom are missing. 

Appendix C: Evaluation of the Vssss integral 

The approximate analytical expression ( flQj) was proposed in Ref. |2| for the two electron 
integral Vssss- By comparing the analytical results with the explicit ab initio ones (dotted 
curve and filled circles in Fig. El respectively) we found a progressively growing deviation 
between them with increasing Z (1.5 eV for lithium, 4.5 eV for fluorine). This deviation can 
be removed by substituting 1.037-Z instead of Z in Eq. (flO]) . where the coefficient 1.037 has 
been extracted from a linear fit (dashed curve in Fig. [6]) of the ab initio results. 

Appendix D: Effect of dynamic correlations 

We performed the CI calculations with the space of the singly excited configurations from 
the CASSCF configurations and examined the semi-internal correlation. The results for the 
single and double core-hole states of C2H2 and CO were summarized in Table IVHI The 
difference between the results of CI and CASSCF provides the effect of the semi-internal 
correlation. The semi-internal correlation has small effect on IPs of the single hole states, 
less than +0.03 eV for Cls -1 and +0.25 eV for Ols -1 . The effect for the one-site double 
core-hole states is significant as +0.24 ~ +0.54 eV, while it is small for the two-site double 
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core-hole states as +0.08 ~ +0.12 eV. 
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TABLE I. Single core hole ionization potentials and their constituent parts (in eV) as calculated 

with the ASCF and CASSCF methods. Experimental values of the ionization potentials, where 

available, taken from Refs. 38-44 are also shown. The cc-pVTZ basis sets were employed. 

Core level, ASCF CASSCF 

Molecule -eg Expt. 

S IP RiS- 1 ) IP RCiS- 1 ) CiS' 1 ) 

LiF Lils 66.407 65.460 0.947 65.334 1.073 0.126 

Fls 710.484 688.187 22.297 688.044 22.440 0.143 

BeO Bels 128.383 124.406 3.977 123.399 4.984 1.007 

Ols 556.694 535.181 21.513 535.075 21.619 0.106 

BF Bis 209.735 202.539 7.196 201.724 8.011 0.815 

Fls 717.810 695.873 21.937 695.915 21.895 -0.042 

CO Cls 309.111 298.256 10.855 296.358 12.753 1.898 296.069 

Ols 562.348 542.801 19.547 542.820 19.528 -0.019 542.543 

N 2 N^s 427.159 409.615 17.544 408.614 18.545 1.001 408.5 

N c ls 432.005 415.373 16.632 412.524 19.481 2.849 412.5 

Ols 563.760 543.046 20.714 542.537 21.223 0.509 542.0 

C0 2 Cls 311.930 300.607 11.323 297.647 14.283 2.960 296.78 

Olg 561.956 19.977 542.870 19.086 -0.891 

541.979 540.6 

Olu 561.954 19.975 542.868 19.086 -0.889 

C 2 H 2 Clg 305.897 13.362 292.202 14.062 0.700 

292.535 291.14, 291.20 

Clu 305.794 13.259 292.111 14.054 0.795 

C 2 H 4 Clg 305.557 13.756 291.344 14.213 0.457 

291.801 290.70, 290.88 

Clu 305.508 13.707 291.297 14.211 0.504 

C 2 H 6 Clg 305.040 13.266 291.147 13.893 0.627 

291.774 290.76 290.74 

Clu 305.023 13.249 291.125 13.898 0.649 

N 2 Nig 426.686 15.444 411.027 15.659 0.215 409.93 

411.242 

Nlu 426.588 15.346 410.932 15.656 0.310 409.82 
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TABLE II. Calculated double core hole ionization potentials and the static correlation energies C 
(in eV) . T and S refer to triplet and singlet couplings of two core holes created on different atomic 
sites, respectively. The cc-pVTZ basis sets were employed. 
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C0 2 Cls- 2 
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TABLE III. Kinetic energies of photoelectrons one would detect in an XTPPS experiment by 
irradiating the CO molecule with an x-ray pulse with photon energies of 1 keV. The sequence of 
the core vacancies in the records reflects the sequence of ionization steps. T and S refer to triplet 
and singlet couplings of two core holes created on different atomic sites, respectively. 

State KE (eV) 

Cls- 1 703.642 

Cls^Cls- 1 631.940 

Ols^Cls" 1 , S 688.077 

Ols^Cls- 1 T 687.620 

Ols- 1 457.180 

Ols^Ols- 1 366.259 

Cls^Ols -1 , S 441.615 

Cls^Ols -1 , T 441.518 
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TABLE IV. Calculated energy differences AE1(S~ 2 ) and AE2(S~ 1 ,S j x ) together with the intra- 
and interatomic generalized relaxation energies RC(S~ 1 ) and IRC(S^ 1 , S^ 1 ) (in eV). 



Molecule 


Energy difference Generalized relaxation 


energy 


LiF 


A£l(Lils- 2 ) 


41.927 ^(Lils- 1 ) 


1.41 




A£^l(Fls" 2 ) 


105.407 ^(Fls- 1 ) 


20.17 




A-E^Lils -1 , Fls- 1 ) 


9.833 LRC*(Lils-\ Fls" 1 ) 


-0.74 


BeO 


A£l(Bels" 2 ) 


51.234 ^(Bels- 1 ) 


5.23 




A£l(01s~ 2 ) 


89.201 ^(Ols" 1 ) 


19.90 




AE2(Bels-\ Ols- 1 ) 


13.327 IRC(Bels-\ Ols" 1 ) 


-2.69 


BF 


A£1(B1 S - 2 ) 


61.875 iJCXBls- 1 ) 


8.38 




AEl(Fls- 2 ) 


103.979 ^(Fls- 1 ) 


20.88 




A J B2(Bls- 1 , Fls" 1 ) 


12.929 LRC*(Bls-\ Fls- 1 ) 


-1.52 


CO 


A£1(C1 S - 2 ) 


71.702 ^(Cls" 1 ) 


11.87 




A£l(01s- 2 ) 


90.921 ^(Ols- 1 ) 


19.03 




A J B2(Cls~ 1 , Ols- 1 ) 


15.565 IRCiCls' 1 , Ols- 1 ) 


-2.80 


N 2 


AE1(N1 S - 2 ) 


79.196 iJCXNls- 1 ) 


17.65 




A£2(Nils-\ N2IS- 1 ) 


13.825 /i2C(Nils _1 , Nsls" 1 ) 


-0.65 


C2H2 


A£l(Cls- 2 ) 


66.653 iiCXCls- 1 ) 


15.06 




A J B2(Cils- 1 , C2IS- 1 ) 


11.015 IRCidla- 1 , C2IS- 1 ) 


0.96 


C2H4 


A£l(Cls- 2 ) 


65.915 iiCXCls- 1 ) 


15.42 




A J B2(Cils- 1 , C2IS- 1 ) 


8.873 IRCXdls- 1 , C2IS- 1 ) 


1.94 


C2H6 


A£l(Cls- 2 ) 


66.555 iiCXCls- 1 ) 


15.06 




AE2(dls-\ C2IS- 1 ) 


6.735 IRCidls' 1 , C2IS- 1 ) 


2.72 


C0 2 


A J Bl(Cls- 2 ) 


69.335 iJCXCls- 1 ) 


13.76 




A£l(01s- 2 ) 


86.171 ^(Ols- 1 ) 


22.70 




A J B2(Cls- 1 , Ols" 1 ) 


10.572 IRCiCls' 1 , Ols- 1 ) 


1.79 




A J B2(01s- 1 , Ols" 1 ) 


8.352 LRC(01s-\ Ols" 1 ) 


-2.19 


N 2 


A£l(N t ls- 2 ) 


76.698 ^(Ntls- 1 ) 


17.78 




A£1(N C 1 S - 2 ) 


77.264 ^(Ncls" 1 ) 


17.47 




A£1(01 S - 2 ) 


88.175 ^(Ols- 1 ) 


20.40 




A£2(N t ls-\ NJs" 1 ) 


11.827 LRC(N t ls-\ NJs" 1 ) 


1.11 




AE2(N t ls- 1 , Ols" 1 ) 


2fo.732 IRCiNtls' 1 , Ols- 1 ) 


0.07 




A J B2(N c ls- 1 , Ols" 1 ) 


11.890 IRCQicla- 1 , Ols" 1 ) 


-6.00 



TABLE V. Single and double core hole ionization potentials of CO and BF calculated with the 
CASSCF method using the relativistic DKH8 Hamiltonian and the differences between them and 
the respective non-relativistic results given in Tables HI and ITT1 All energies are in eV. T and S refer 
to triplet and singlet couplings of two core holes created on different atomic sites, respectively. 



Molecule 




State 




Ionization Potential Difference 


CO 




Cls -1 




one AAP. 
ZyO.440 


U.Uoo 






Ols -1 




D4O.100 


U.O40 






Cls" 2 




004.001 


U.Z4o 






01s- 2 




1177.408 


0.847 




Cls 


^Ols -1 . 


S 


855.176 


0.433 




Clfl- 


^Ols- 1 , 


T 


855.632 


0.432 


BF 




Bis" 1 




201.758 


0.034 






Fls" 1 




696.500 


0.585 






Bis" 2 




465.425 


0.102 






Fls" 2 




1497.177 


1.368 




Bls 


-^Fls" 1 , 


S 


911.187 


0.619 




Bis 


^Fla" 1 , 


T 


911.297 


0.619 
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TABLE VI. Basis set dependence of the CASSCF single and double core hole ionization potentials. 
All energies are in eV. T and S refer to triplet and singlet couplings of two core holes created on 
different atomic sites, respectively. 



Molecule 




State 




cc-pVDZ 


cc-pVTZ 


cc-pCVTZ 


cc-pVQZ 


C2H2 








293.262 


292.202 


292.111 


291.728 






Clu- 1 




293.172 


292.111 


292.020 


291.634 






Cls- 2 




658.167 


650.228 


649.412 


649.711 




CiLs 


- 1 C 2 ls 


-\s 


598.350 


594.590 


594.337 


593.908 




Cils 


~ 1 C 2 ls 




598.853 


595.197 


594.971 


594.592 


CO 




Cls- 1 




298.062 


296.358 


296.239 


296.229 






Ols- 1 




544.881 


542.820 


542.611 


542.559 






Cls- 2 




672.908 


664.418 


663.663 


663.632 






Ols- 2 




1184.846 


1176.561 


1175.469 


1175.373 




Cls 


^Ols- 


\s 


859.146 


854.743 


854.452 


854.425 




Cls 


-iQls- 




859.491 


855.200 


854.931 


854.905 
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TABLE VII. Effect of dynamic correlations examined by the CI calculation. All energies are in eV. 



Molecule 




State 




CASSCF 


CI 


C2H2 




Clg- 1 




292.202 


292.216 






Clu- 1 




292.111 


292.127 










650.228 


650.586 




Cila 


- 1 C 2 l5- 1 


, s 


594.590 


594.728 




Cils 


- 1 c 2 i s - 1 


,T 


595.197 


595.319 


CO 




Cls" 1 




296.358 


296.386 






01s- 1 




542.820 


543.070 






Cls" 2 




664.418 


664.658 






Ols- 2 




1176.561 


1177.096 




Cls 


^Ola" 1 , 


S 


854.743 


854.833 




Cls 


^Ols- 1 , 


T 


855.200 


855.276 
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FIG. 1. (Color online). Schematic picture of x-ray two-photon photoelectron spectroscopy 
(XTPPS) and x-ray two-photon-induced Auger electron spectroscopy (XTPAES). See text for 
explanations. In this picture it is assumed that the second photon is absorbed before Auger decay 
takes place and that one core hole decays much faster than the other. In reality all processes 
overlap. 

FIG. 2. (a) The energy difference AE1(S- 2 ) = Z>JP(5-\ as a function 

of the atomic number Z of the atom with the core orbital S; (b) The generalized relaxation energy 
RC(S^ 1 ) calculated by means of Eq. (fT5j) as a function of Z. 

FIG. 3. The ratio n = RC(S^ 1 , S^ 1 )/ RC(S^ 1 ) as a function of the atomic number Z of the atom 
with the core orbital S. 

FIG. 4. (a) The energy difference AE2(Sr 1 ,S~ 1 ) = DIP(Sr\ Sr 1 ) - [IPiSr 1 ) + IP^Sr 1 )} as 
a function of the inverse distance r between the atoms with the core orbitals Si and Sj; (b) The 
interatomic generalized relaxation energy IRCiS^^Sr 1 ) calculated by means of Eq. (|16j) as a 
function of r. 

FIG. 5. Differences between the valence electron density distributions of the ground and double 
core hole states: (a) Cils" 2 of C 2 H 4 ; (b) Cil S - 1 C 2 ls- 1 (singlet) of C 2 H 4 ; (c) 01s~ 2 of CO 

FIG. 6. The two-electron integral Vssss as a function of the atomic number Z of the atom with the 
core orbital S. The results extracted from the ab initio calculations (filled circles) are compared 
with those calculated by means of Eq. (jlOp (dotted curve). A linear fit of the ab initio results is 
also shown (dashed curve). 

FIGURE CAPTIONS 
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